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CN ■ We perform a detailed investigation of total lifetimes for the doubly heavy baryons 

Hqq/ , 0,qqi in the framework of operator product expansion over the inverse heavy quark 

mass, whereas, to estimate matrix elements of operators obtained in OPE, approxima- 

Q\ • tions of nonrelativistic QCD are used. 
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At present a number of powerful techniques based on Operator Product Expansion (OPE) 
and effective field theories have been developed. These tools allow one consistently to include 
into consideration various nonperturbative contributions, written in terms of a few number of 
universal quantities. The coefficients (Wilson coefficients) in front of these operators are gen- 
erally expanded in series over the QCD coupling constant, inverse heavy quark mass and/or 
relative velocity of heavy quarks inside the hadron. The accuracy, obtained in such calcula- 
tions, can be systematically improved, and it is limited only by the convergence properties of 
the mentioned series. The described approach have been already widely used for making the 
precise predictions in the heavy quark sector of Standard Model (SM), such as decays, dis- 
tributions and partial width asymmetries involving the CP violation^ for the heavy hadrons. 
The sensitivity of Wilson coefficients to virtual corrections caused by some higher-scale inter- 
actions makes this approach to be invaluable in searching for a "new" physics at forthcoming 
experiments. 

The approach under discussion has been successfully used in the description of weak decays 
of the hadrons containing a single heavy quark, as carried out in the framework of Heavy Quark 
Effective Theory (HQET) ||, in the annihilation and radiative decays of heavy quarkonia QQ, 
where one used the framework of non-relativistic QCD (NRQCD) H, and in the weak decays 

1 For review see PJl. 
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of long-lived heavy quarkonium with mixed flavours |4|] |. The experimental data on the 
weak decays of heavy hadrons can be used for the determination of basic properties of weak 
interactions at a fundamental level, in particular, for the extraction of CKM matrix elements. 
The same approach is also valid for the baryons containing two heavy quarks. 

In addition to the information extracted from the analysis of hadrons with a single heavy 
flavor, the baryons with two heavy quarks, QQ'q, provide a way to explore the nonspectator 
effects, where their importance is increased. Here we would like to note, that in the case of 
systems with two heavy quarks, the hypothesis on the quark-hadron duality is more justified, 
and, so, the results of OPE-based approach turn out to be more reliable. For these baryons we 
can apply a method, based on the combined HQET-NRQCD techniques |2|, |], d, if we use the 
quark-diquark picture for the bound states. The expansion in the inverse heavy quark mass for 
the heavy diquark QQ' is a straightforward generalization of these techniques in the mesonic 
decays of B c [|], with the difference that, instead of the color-singlet systems, we deal with 
the color-anti-triplet ones, with the appropriate account for the interaction with the light 
quark. First estimates of the lifetimes for the doubly heavy baryons and S£ c were recently 
performed in [^, ||. Using the same approach, but different values of parameters^ repetition 
of our results for the case of doubly charmed baryons was done in [|K|. The spectroscopic 
characteristics of baryons with two heavy quarks and the mechanisms of their production in 
different interactions were discussed in refs. ffl], [L2| [L3|, [14| and |l5j , respectively. 



In this paper, we present the calculation of lifetimes for the doubly heavy baryons as well as 
reconsider the previous estimates with a use of slightly different set of parameters adjusted in 
the consideration of lifetime data for the observed heavy hadrons and improved spectroscopic 
inputs. As we made in the description of inclusive decays of the and S£ c -baryons, we 
follow the papers j|, |TB| , where all necessary generalizations to the case of hadrons with two 



heavy quarks and other corrections are discussed. We note, that in the leading order of 
OPE expansion, the inclusive widths are determined by the mechanism of spectator decays 
involving free quarks, wherein the corrections in the perturbative QCD are taken into account. 
The introduction of subleading terms in the expansion over the inverse heavy quark massesf] 
allows one to take into account the corrections due to the quark confinement inside the hadron. 
Here, an essential role is played by both the motion of heavy quark inside the hadron and 
chromomagnetic interactions of quarks. The important ingredient of such corrections in the 
baryons with two heavy quarks is the presence of a compact heavy diquark, which implies that 
the square of heavy quark momentum is enhanced in comparison with the corresponding value 
for the hadrons with a single heavy quark. The next characteristic feature of baryons with 
two heavy quarks is the significant numerical impact on the lifetimes by the quark contents 
of hadrons, since in the third order over the inverse heavy quark mass, 1/mg, the four-quark 
correlations in the total width are enforced in the effective lagrangian due to the two-particle 
phase space of intermediate states (see the discussion in fll6||). In this situation, we have to add 
the effects of Pauli interference between the products of heavy quark decays and the quarks in 
the initial state as well as the weak scattering involving the quarks composing the hadron. Due 



2 The first experimental observation of the B c -meson was recently reported by the CDF Collaboration ||; 
see ref. || for a theoretical review of £? c -meson physics before the observation. 

3 See comments on the difference in the numerical values of lifetimes of doubly charmed baryons in [pj. 

4 It was shown in |L7| that the first order 1/mQ-correction is absent, and the corrections begin with the 
1/mQ-terms. 
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to such terms we introduce the corrections depending on spectators and involving the masses 
of light and strange quarks in the framework of non-relativistic models with the constituent 
quarks, because they determine the effective physical phase spaces, strongly deviating from 
the naive estimates in th decays of charmed quarks. We take into account the corrections 
to the effective weak lagrangian due to the evolution of Wilson coefficients from the scale of 
the order of heavy quark mass to the energy, characterizing the binding of quarks inside the 
hadron |16| . 

The paper is organized as follows. In agreement with the general picture given above, we 
describe the scheme for the construction of OPE for the total width of baryons containing 
two heavy quarks with account of corrections to the spectator widths in Section 2. The 
procedure for the estimation of non-perturbative matrix elements of operators in the doubly 
heavy baryons is considered in Section 3 in terms of non-relativistic heavy quarks. Section 4 
is devoted to the numerical evaluation and discussion of parameter dependence of lifetimes of 
doubly heavy baryons. We conclude in section 5 by summarizing our results. 



2 Description of the method 

In this section we describe the approach used for the calculation of total lifetimes for the doubly 
heavy baryons, originally formulated in |7|, ||, together with some new formulae, required for 
the evaluation of nonspectator effects in the decays of otheiQ baryons in the family of doubly 
heavy baryons, not considered previously. 

The optical theorem along with the hypothesis of integral quark-hadron duality, leads us to 
a relation between the total decay width of heavy quark and the imaginary part of its forward 
scattering amplitude. This relationship, applied to the S^L-baryon total decay width r„(») , 

"QQ' 

can be written down as: 

1 
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where the Hqq, state in Eq. ([[]) has the ordinary relativistic normalization, (Hqq,|5qq,) 
2EV, and the transition operator T is determined by the expression 



T = %m <fx {TH eff (x)H eff (0)}, 



(2) 



where H e ff is the standard effective hamiltonian, describing the low energy interactions of 
initial quarks with the decays products, so that 
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where f is the number of flavors, {a, /3, 7, 5} run over the color indeces. 

Under an assumption, that the energy release in the heavy quark decay is large, we can 
perform the operator product expansion for the transition operator T in Eq.(|I]). In this 
way we obtain series of local operators with increasing dimensions over the energy scale, 
wherein the contributions to T„(,) are suppressed by the increasing inverse powers of the 

"QQ' 

heavy quark masses. This formalism has already been applied to calculate the total decay 
rates for the hadrons, containing a single heavy quark JlTj] (for the most early work, having 
used similar methods, see also |TjJ [L8|]) and hadrons, containing two heavy quarks [0, ||. As 



was already pointed in j7|, the expansion, applied here, is simultaneously in the powers of 
both inverse heavy quark masses and the relative velocity of heavy quarks inside the hadron. 
Thus, this fact shows the difference between the description for the doubly heavv baryons and 
the consideration of both the heavy- light mesons (the expansion in powers of QCD } and the 

heavy-heavy mesons |4| (the expansion in powers of relative velocity of heavy quarks inside 
the hadron, where one can apply the scaling rules of nonrelativistic QCD fl9|j ). 
The operator product expansion explored has the form: 

2 1 1 

T = V C 1 (^)Q i Q i + —C 2 {^ga, v G^Q* + —0(1) (4) 

The leading contribution in Eq.(f|) is determined by the operators Q l Q t , corresponding 
to the spectator decay of Q l -quarks. The use of motion equation for the heavy quark fields 
allows one to eliminate some redundant operators, so that no operators of dimension four 
contribute. There is a single operator of dimension five, Q 1 G q = Q l go '^G^ 1 'Q* As we will show 
below, significant contributions come from the operators of dimension six Q 2 Q i 2q — Q l TqqT' Q l , 
representing the effects of Pauli interference and weak scattering for doubly heavy baryons. 
Furthermore, there are also other operators of dimension six = Q l a flu , jiD tJ- G ul Q t and 

Qe2Qi — Q 1 D ^G^T V Q % , which are suppressed in comparison with Q 2 Qi 2g [0]- 111 what follows, 
we do not calculate the corresponding coefficient functions for the latter two operators, so that 
the expansion is certainly complete up to the second order of — , only. 

Further, the different contributions to OPE are given by the following: 

?n+ = 2 ^35c + T^p 1, 

be 1 ' 

T~? bc = %5b + T 35c + T^pj + T^ ws , 
%- = 2%5b + Tq%, 
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where the 35- labelled terms account for the operators of dimension three 3 gi and five Ogqm 
the 6-marked terms correspond to the effects of Pauli interference and weak scattering. The 
explicit formulae for these contributions have the following form: 



T 35b = T Kspec bb - ^[2P cl + P CTl + K ob (P cl + P ccl ) + K 2b (P c2 + P cc2 )}0, 



m 



Gb, 



(5) 



where 
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2C 2 , K 2 q = 2{C+ — C£), and Tq spec denotes the spectator width (see 
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where y — and r = m 2 /m 2 . The functions P s i(P s2 ) can be obtained from P c i(P C 2) by the 

substitution y = m 2 /m 2 . In the 6-quark decays, we neglect the value m 2 /m 2 and suppose 
m s = 0. 

The calculation of both the Pauli interference effect for the products of heavy quark decays 
with the quarks in the initial state and the weak scattering of quarks, composing the hadron, 
depends on the quark contents of baryons and results in: 



(i) 
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^6,pi — 1~p I ud + ^2 Tpi, Ul i + ?~pi,sc + Tp I & + ^ ^pim + Tp I sS (19) 
i i 

Tq%s = T W s,bc + Tws,cs (20) 

^!ws = ^wsm (21) 

^pj = 2Tp b Ida (22) 

= 27^ (23) 

G 2 F \V cb \ 2 2 m c 2 
-m b (l 



([(^-^ " - 75)&0(97 Q (1 - 75)9) + 

( (1 ~/" )2 " (1 ~ 3 ^ )3 )(^75^)(c J 7 a (l - 75)9)] (24) 
[(C+ - C_) 2 + ^(1 - A;5)(5C 2 + d + 6C_C + )] + 

[(^^ " { ^—^){hl a {l - 76)6 J -)(S,-7 a (l - 7 5 )Q) + 

( (1 ~/" )2 - (1 ~3^ )3 )(^75^)(c j7 a (l - 75 ) Ci )]^(5C 2 + C! + 6C_C + )), 

^IjKbP 2/1 m ^2 

\^-^~ " - 75)fe J )(c j 7 a (l - 75 )Q) + (25) 

(^-^ - S-^)(5 i 7a7fe6i)(S,-7 a (l -7b)ci)], 
^'^^-m^l - — ) 2 



47r m b ' 

([(2^^ " ^^W(l - 75)^)K'7 a (l - 75K) + 
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[(^^ - 2-^)(^7a(l - 7 5 )^)K-7 Q (1 - 75)*) + 

( (1 ~/" )2 - (1 V" )3 )(5<7a756j)(Jj7a(1 " 76)*)]*'(5CJ + C 2 - 6C_C + )), 

l07T 171b 

([(1 - z_)(& i7a (l - 7s)^)(^7 a (l - 75)^) + 

2 -{l + 2z_)(fo i7 a75&i)(^7 Q (l - 75)^)] (27) 
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(c i7 a(l - 75)Q)(^7 a (l - 75)^) + (33) 
kHcl - C 2 _)(c lla (l ~ 75b)K-7 a (l - 75)^)], 
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where the following notation has been used: 
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In the evolution of coefficients C+ and C_, we have taken into account the threshold effects, 
connected to the heavy quark masses. 

In expressions fl5|) and @, the scale /1 has been taken approximately equal to m c . In 
the Pauli interference term, we suggest that the scale can be determined on the basis of the 
agreement of the experimentally known difference between the lifetimes of A c , and S° 
with the theoretical predictions in the framework described abovef]. In any case, the choice of 
the normalization scale leads to uncertainties in the final results. Theoretical accuracy can be 
improved by the calculation of next-order corrections in the powers of QCD coupling constant. 

The coefficients of leading terms, represented by operators bb and cc, are equivalent to 
the widths fot the decays of free quarks and are known in the next-to-leading logarithmic 
approximation of QCD |23|, ^4| |25|, [27j, including the mass corrections in the final state 

6 A more extended description is presented in Q| . 
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with the charmed quark and r-lepton [27j in the decays of 6-quark and with the strange quark 
mass for the decays of c-quark. In the numerical estimates, we include these corrections and 
mass effects, but we neglect the decay modes suppressed by the Cabibbo angle, and also the 
strange quark mass effects in b decays. 

The expressions for the contribution of operator Ym=i ^GQ i are known in the leading loga- 
rithmic approximation. The expressions for the contributions of operators with the dimension 
6 have been calculated by us with account of masses in the final states together with the 
logarithmic renormalization of the effective lagrangian for the non-relativistic heavy quarks at 
energies less than the heavy quark masses. 

Thus, the calculation of lifetimes for the baryons Sqq, is reduced to the problem of evalu- 
ating the matrix elements of operators, which is the subject of next section. 



3 Matrix elements in NRQCD approximation. 

By using the equations of motion, the matrix element of operator QiQi can be expanded in 
series over the powers of l/rngj: 

l*> \ -1 ^Q Q m^D) 2 -(^G)]Q^ QQI ) norm l 

zln Qj " L Qi 

Thus, we have to estimate the matrix elements of operators from the following list only: 

Qi(iDfQi, (~)Q j aGQ j , Q^ a (l - 75 )Q% a (l - 7s)?, 

a'',,'.-/>V''(l - l 5 )q, Q j lal 5 Q j Q k l a {l - l 5 )Q k , (38) 
Q J 7a(l-75)Q J QV(l-75)Q fe . 

The meaning of each term in the above list, was already discussed by us in the previous papers 
on the decays of doubly heavy baryons |7], H|, so we omit it here. 

Further, employing the NRQCD expansion of operators QQ and Qgcr^G^Q, we have 

QQ 



Qg^GTQ 

Here the factorization at scale \i (tuq > fi > tjiqVq) is supposed. We have omitted the term of 
Wqa{gE x D)^/q, corresponding to the spin-orbital interactions, which are not essential for 
the basic state of baryons under consideration. The field \1/q has the standard non-relativistic 
normalization. 

Now we would like to make some comments on the difference between the descriptions 
of interactions of the heavy quark with the light and heavy heavy ones. As well known, in 
the doubly heavy subsystem there is an additional parameter which is the relative velocity of 



l -T^ Q 9° B *Q ~ ^q( D 9E)^q + ••• (39) 



-2tf Q gtrB* Q - —^ Q (DgE)^ Q + ... (40) 

m Q 
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quarks. It introduces the energy scale equal to ttiqV. Therefore, the Darwin term {DE) in 
the heavy subsystem stands in the same order of inverse heavy quark mass in comparison with 
the chromomagnetic term (<tB) (they have the same power in the velocity v) . This statement 
becomes evident if we apply the scaling rules of NRQCD 



3 — 

^Q^(m Q v Q )^, \D\~m Q v Q , gE ~ m 2 Q v^, gB ~ m 2 Q v Q , g ~ v%. 

For the interaction of heavy quark with the light one, there is no such small velocity parameter, 
so that the Darwin term is suppressed by the additional factor of k/rriQ ~ Kqcd/^q- 

Further, the phenomenological experience with the potential quark models shows, that the 
kinetic energy of quarks practically does not depend on the quark contents of system, and it 
is determined by the color structure of state. So, we suppose that the kinetic energy is equal 
to T = rridV J/2 + mivf/2 for the quark-diquark system, and it is T/2 = m^/2 + m c v 2 /2 in 
the di quark (the color factor of 1/2). Then 

(Zl Q ,\tf Q (iDy* Q \Z« QQI ) _ 2 ^ 2m q T | m Q ,T 



2M^ QQ m 2 Q y (m, + m Q > + m Q )(m Q > + m Q ) m Q {m Q + m Q > 

(Zq Q H^ D ) 2 *Q'\ E QqK _ v 2^_ 1rn q T | m Q T 

2M s > QQ m 2 QI Q ' (m q + m Q > +m Q )(m Q/ +m Q ) m Q >{m Q + m Q /)' 

where the diquark terms dominate certainly. Applying the quark-diquark approximation and 
relating the matrix element of chromomagnetic interaction of diquark with the light quark to 
the mass difference between the exited and ground states M=«* — M=o , we have 

° QQ' QQ' 



2 + ... 



(H<L|cc|S<L) 1 1 2 lM 3 o. -M S o 5g 2 |T/ 
2M^o 2 c 3 m c 18m? 1 

'cc L ' C 

» 1 - 0.073 - 0.025 - 0.009 + ... (43) 

RJsW _ 1 2 _ _ 5^ 2 a 

2Msjo 2 c 3 m c 18m3 1 1 )l 

w 1 - 0.078 - 0.025 - 0.009 + . . . (44) 

""""be L ^ 

« 1 - 0.098 + 0.006 - 0.010... (45) 

^^mp - i -^ + 3^™-^ 2 i* i (°)i 2+ --- 

« 1 - 0.099 + 0.006 - 0.010... (46) 

\ lih\ —///// J — t; "ftfc _ | 'I' (0)| + 



2M H - 6 2 3 m b 18m£ 



» 1 - 0.021 - 0.003 - 0.002 + ... (47) 

WW . _ 1 ,2 _ 1 Mngg - M n o fc _ 5ff 2 

2M fi o 6 2 fe 3 m b I8m 3 b l { >l + 

» 1 - 0.021 - 0.003 - 0.002 + ... (48) 
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The numerical values of parameters used in the calculations above are given in the next section. 
Our presentation here is less detailed than in previous papers |7|, |8|]. However, we hope, that 
the interested reader can find there all needed details. 

Analogous expressions may be obtained for the matrix elements of operator Qgcr^G^Q 



2M=o m 2 r 
(~i c \cga, u G^c\ZZ c ) 
2M=« m 2 r 

^bc c 

(~l b \bga, u G^b\^) 
2M=o m 2 h 

^bb 

{Vt QQ ,\cg(j^G^c\VLQQ,) 



4 (M=o, - M=o 
3 m r 



3m b mt 



I* 



9m? 



-0.124, 



3m;? 



^ d (0)| 2 « 0.005, 



4 (M 5 q 
3 m?, 



7^ 



•qj c 



-0.189, 



(49) 
(50) 

(51) 
(52) 



The permutations of quark masses lead to the required expressions for the operators of bb and 
bga^G^b. 

For the four quark operators, determining the Pauli interference and the weak scattering, 
we use the estimates in the framework of non-relativistic potential model [7|, ||: 



(^1(^(1-75)0)^(1 


-7s)g)|2* c ) 


= 12(m c + m g ) 


■ |# dZ (0)| 2 , 


(53) 




-75)?)|SD 


= 8(m c + m g )- 


|^(0)| 2 , 


(54) 


(ft cc |(c 7M (l-75)c)(s7 M (l 


- 7 5 )s)|fi cc ) 


= 12(m c + m s ) 


• |^(0)| 2 , 


(55) 


(fi cc |(c7 M 7 5 c)(s7^(l 


- 7 5 )s)|n cc ) 


= 8(m c + m s )- 


|^(0)| 2 , 


(56) 


(~l\(h,(l- l5 )b)(cr(l 


-7b)c)|S^) 


= 8(m b + m c )- 


|^ d (0)| 2 , 


(57) 




-75)c)|^ c ) 


= 6(m b + m c )- 


|^ d (0)| 2 , 


(58) 


(El\(h,(l-l 5 )b)(qY(l 


-75)g)|^ c ) 


= 2(m b + mi)- 


|^(0)| 2 , 


(59) 


(KJiGhsrsb) (gy(i 


-7s)g)l^ c ) 


= 0, 




(60) 


(S^|(c7 At (l-75)c)(g 7 ' 1 (l 


-7s)g)|2fe c ) 


= 2(m c + mi)- 


|^(0)| 2 , 


(61) 


(HLI(c7 M 75c)(g 7 M (l 


-7 5 )g)|HL) 


= 0, 




(62) 


(Q bc \(bj,(l- 75)6) (cr^l 
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= 8(m b + m c )- 


|^(0)| 2 , 


(63) 


(0 6c |(67 M 7 5 &)(c 7 ^(l 


- 75 )c)|0 6c ) 


= 6(m b + m c )- 


|^ d (0)| 2 , 


(64) 


(^1(67^(1 -7 5 )&)(S7 M (1 


- 75)s) \ttbc) 


= 2(m b + m s )- 


|^(0)| 2 , 


(65) 


(^ C |(&7/,75&)(^(1 


- 75) s) n 6c ) 


= 0, 




(66) 


(Q 6c |(c7 M (l-7 5 )c)(s7^(l 


- 75)s)|0 6c ) 


= 2(m c + m s )- 


|^(0)| 2 , 


(67) 


(fi bc |(c7 M 7 5 c)(s7^(l 


- 7 5 )s) |0 6c ) 


= 0, 




(68) 


(^ 6 |(6 TM (l-75)fc)(g7 M (l 


-7B)g)|si> 


= 12(m b + m q ) 


■ |^'(0)| 2 , 


(69) 


(Ht b |(&7 M 75&)(g7 M (l 


-75)g)|2fe 6 ) 


= 8(m b + m q )- 


|^(0)| 2 . 


(70) 



The color structure of wave functions leads to the relations 

{^ QQ ,mTM{q k ^{i - 7fe)ft)|s^> = -(^ QQ/ \(QT,Q)(mi - 75)g)|^ Q ,>, 

where T M is an arbitrary spinor matrix. 
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4 Numerical estimates 

Performing the numerical calculations of lifetimes for the doubly heavy baryons, we have used 
the following set of parameters: 

m s = 0.2 GeV m, = 0. GeV m* = 0.45 GeV m{ = 0.3 GeV 

\V CS \ = 0.9745 \V bc \ = 0.04 T = 0.4 GeV (71) 
m c = 1.55 GeV m b = 5.05 GeV (72) 

The numerical values of diquark wavefunctions at the origin for baryons under consideration 
are collected in Table [I]. The masses of doubly heavy baryons may be found in Table ||. 





-^++ 

' — 'cc 


' — 'cc 


^tc 


^bc 


^0 
"6c 


0° 

"fee 


^0 

"fefe 


"66 


^66 


tf d (0), GeV-2 


0.150 


0.150 


0.150 


0.205 


0.205 


0.205 


0.380 


0.380 


0.380 



Table 1: The values of diquark wavefunctions for the doubly heavy baryons at the origin. 





-++ 

' — 'cc 


' — 'cc 


cc 


"6c 


^0 

"6c 




^0 
"66 


"66 


^6~6 


M, GeV 


3.478 


3.478 


3.578 


6.82 


6.82 


6.92 


10.093 


10.093 


10.193 


M*, GeV 


3.61 


3.61 


3.71 








10.193 


10.193 


10.293 



Table 2: The masses of doubly heavy baryons M, and M* stands for the mass of the baryon 
with lowest excited state of light quark-diquark system. 

The wavefunctions as well as masses for the considered baryons are taken from |i~2] , where 
their estimates in the non-relativistic model with the Buchmiiller-Tye potential were done. 
The 6-quark mass is obtained from the requirement, that for any given value of c-quark mass 
the theoretically computed I^-meson lifetime equals to experimentally measured value. This 
matching condition leads to the following approximate relation 



3.5 GeV. 



(73) 



The c-quark mass is determined from the analogous matching procedure for the 5 c -meson 
lifetime 0] . The m*-values in Eq. ( ffTD represent the constituent masses for the corresponding 
light quarks, used by us in estimations of hadronic matrix elements.^] For the value of light 
quark-diquark function at the origin we assume 



2p D M D k~ 
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(74) 



where fn = 170 MeV. This expression obtained by performing the steps similar to [28, 2£ 



for the derivation of hyperfine splitting in the light quark-diquark system. The factor k~ 



'See M for details. 
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1 — 'cc 


1 — cc 


n+ 




3.104 


3.104 


3.104 


PI, ps" 1 


-0.874 




0.621 


WS, ps" 1 




1.776 




r, ps 


0.45 


0.20 


0.27 



Table 3: The lifetimes of doubly charmed baryons together with the relative spectator and 
nonspectator contributions to the total widths. 



accounts for the low energy logarithmic renormalization of fp constant. We have used this 
equation for all doubly heavy baryons, considered in this paper. Even though we use this 
relation to compute central values of lifetimes, the precise values of wavefunction parameters 
are under question, so in the presented results we have allowed for variations. 

The renormalization scale fi is chosen in the following way: fii = rrib and fi2 = itl c in the 
estimates of Wilson coefficients C+(/x) and C-(p) for the effective four-fermion weak lagrangian 
at low energies with the b and c-quarks, correspondingly. For nonspectator effects, which are 
the Pauli interference and weak scattering of valence quarks, the renormalization scale \x is 
obtained from the fit of theoretical predictions for the lifetime differences of baryons A c , 
and S° over the experimental data. 

In Table § we present the results of calculations for the doubly charmed baryons. Together 
with the total lifetimes of these baryons we have shown the relative spectator and nonspectator 
contributions. From this Table we see the importance of nonspectator effects, producing huge 
differences in the values of lifetimes. The analogous results for other doubly heavy baryons 
can be found in Tables [| and |. 

A small comment concerns with the corrections to the spectator decays of heavy quarks, 
caused by the motion of heavy quarks inside the hadron and interactions with the light degrees 
of freedom. The corrections due to the quark-gluon operators of dimension 5 are numerically 



small ]TJ| . The most important terms come from the kinetic energy of heavy quarks. 



In Figs. we have shown the dependence of baryons lifetimes from the values of light 
quark-diquark wavefunctions at the origin. We see quite a different behaviour withthe increase 
of | ^^(0) [-parameter. 





w+ 

^bc 


wo 

"6c 


0° 

"fee 


b -»• c, ps 1 


0.632 


0.632 


0.631 


Y^c -> s, ps" 1 


1.511 


1.511 


1.509 


PI, ps" 1 


0.807 


0.855 


0.979 


WS, ps" 1 


0.653 


0.795 


1.713 


r, ps 


0.28 


0.26 


0.21 



Table 4: The lifetimes of (fecg)-baryons together with the relative spectator and nonspectator 
contributions to the total widths. 
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^0 
"66 


1=1 66 


^66 


X> -> c, ps" 1 


1.254 


1.254 


1.254 


PI, ps" 1 




-0.0130 


-0.0100 


WS, ps" 1 


0.0189 






r, ps 


0.79 


0.80 


0.80 



Table 5: The lifetimes of (66g)-baryons together with the relative spectator and nonspectator 
contributions to the total widths. 

Here, we would like to note, that in this paper we do not give a detail discussion of 
nonspectator effects on the total lifetimes and semileptonic branching ratios of doubly heavy 
baryons and promise to fill this gap in one of our subsequent papers f30|j . 

Finally, concerning the uncertainties of the presented estimates, we note that they are 
mainly determined by the following: 

1) The c-quark mass is poorly known, but constrained by the fits to the experimental data, 
discussed above, can lead to the uncertainty ~ 15% in the case of doubly charmed baryons 
and t£ ~ 10% for the case of bcq - baryons. 

2) The uncertainties in the values of diquark and light quark-diquark wavefunctions lead 
to jr ~ 30% in the case of doubly charmed baryons and ^ ~ 15% for the bcq - baryons. 

Thus, the estimated uncertainty in predictions for the lifetimes of doubly heavy baryons 
is close to 25% in the case of (bcq) - baryons, of order of 45% in the case of doubly charmed 
baryons and less then 5% in the case of ipbq) - baryons. 

T s ++, pS 



. 65 - 

. 6 : 
. 55 : 

. 5 : 
. 45 : 

0.4^ 

0.002 0.0025 0.003 0.0035 0.004 0.0045 0.005' 

Figure 1: The dependence of H^-baryon lifetime on the value of wavefunction of light quark- 
diquark system at the origin |\l/ di (0)|. 
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^(0)| 2 , GeV 3 



0.002 0.0025 0.003 0.0035 0.004 0.0045 0.005 

Figure 2: The dependence of H+-baryon lifetime on the value of wavefunction of light quark- 
diquark system at the origin |\l/ d '(0)|. 




. 25 



0.24- 



0.002 0.0025 0.003 0.0035 0.004 0.0045 0.005 



|^(0) | 2 , GeV 3 



Figure 3: The dependence of f2 cc -baryon lifetime on the value of wavefunction of light quark- 
diquark system at the origin |\l/ d '(0)|. 
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T s + , pS 

be 



.28 



0.279 



0.278 



.277 



0.002 0.0025 0.003 0.0035 0.004 0.0045 0.005 



|^(0) | 2 , GeV 3 



Figure 4: The dependence of S^,-baryon lifetime on the value of wavefunction of light quark- 
diquark system at the origin |^ d '(0)|. 



T=o , PS 

be 



0.266 



.264 



.262 - 



0.26- 



0.258 - 



. 256 



0.002 0.0025 0.003 0.0035 0.004 0.0045 0.005 



|^'(0)| 2 , GeV 3 



Figure 5: The dependence of S° c -baryon lifetime on the value of wavefunction of light quark- 
diquark system at the origin |^ di (0)|. 
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0.19 



0.18 



0.002 0.0025 0.003 0.0035 0.004 0.0045 0.005 



|^(0) | 2 , GeV 3 



Figure 6: The dependence of f2;, c -baryon lifetime on the value of wavefunction of light quark- 
diquark system at theorigin |^ d '(0)|. 




diquark system at the origin \fy dl (0) 



17 




diquark system at the origin |\l/ dZ (0)|. 



T n bb , PS 




0.002 0.0025 0.003 0.0035 0.004 0.0045 0.005 



|^(0) | 2 , GeV 3 



Figure 9: The dependence of f^-baryon lifetime on the value of wavefunction of light quark- 
diquark system at the origin |^ d '(0)|. 
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5 Conclusion 



In the present paper we have performed a detail investigation and numerical estimates for 
the lifetimes of doubly heavy baryons. The used approach is based on OPE expansion of 
total widths for the corresponding hadrons, and it is combined with the formalism of effective 
fields theories developed previously In this way, we have accounted for the both perturbative 
QCD and mass corrections to the Wilson coefficients of operators. The nonspectator effects, 
presented by Pauli interference and weak scattering, and their influence on the total lifetimes 
are considered. The obtained results show the significant role played by them in the description 
of lifetimes of doubly heavy baryons. 
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tation of results. 
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